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Intisari 

Ultrasonografi (USG) cukup banyak digunakan dalam imaging diagnostic , oleh karena itu kontrol kualitas dari 
citra USG menjadi hal yang penting, salah satu parameternya terkait resolusi spasial. Resolusi spasial dapat di- 
analisis dengan menghitung jarak antar objek terdekat pada arah tangensial citra USG. Citra yang dihasilkan dari 
transduser dengan jenis curved array akan menghasilkan citra yang berbentuk seperti potongan lingkaran dan 
koordinat objek yang dihasilkan dari transduser dengan jenis ini tidak dalam sistem koordinat kartesian karena 
menyesuaikan dengan permukaan dari transduser, sehingga perlu dilakukan transformasi dengan pendekatan 
persamaan lingkaran untuk dapat mengukur jarak antar objek pada arah tangensial. Penelitian dilakukan meng¬ 
gunakan pesawat USG 3D merk Mindray model: DP-10 dengan transduser model: 35C50EB, serta imstrumen 
phantom yang berasal dari campuran silicon rubber 99 ml dan katalis 1 ml. Akuisisi data diambil pada frekuensi 
4,5 MHz, gain sebesar 168 dB, dan depth sebesar 3,3 cm. Setelah dilakukan transformasi menggunakan pen¬ 
dekatan persamaan lingkaran, didapati hasil pengukuran objek pada daerah A sebesar (0,85 ±0,01) cm, daerah 
B sebesar (0,62 ± 0,07) cm, dan daerah C sebesar (0,48 ± 0,06) cm. Nilai-nilai yang didapatkan tersebut masih 
memenuhi standar yang ditetapkan American Association of Physicist in Medicine (AAPM). 

Abstract 

Ultrasound is widely used in diagnostic imaging, therefore quality control (QC) of ultrasound images is impor¬ 
tant. One of the QC parameters is a spatial resolution that can be analyzed by calculating the distance between 
two nearby objects in the tangential direction of radiation. The image that is produced by the curved array trans¬ 
ducer will produce a circular shaped image and the coordinates of the object arent in the cartesian coordinate 
system because it adapts to the surface of the transducer. The image must be transformed with a circle equation 
approach to calculate the distance between two nearby objects. This study was using a Mindray 3D ultrasound 
model: DP-10 with transducer model: 35C50EB, and phantom instruments made from a mixture of 99 ml sili¬ 
con rubber and 1 ml catalyst. The acquisition data used a fixed frequency of 4.5 MHz, fixed-gain of 168 dB, and 
fixed depth of 3,3 cm. The result found that the measurement of objects in A area was (0.85 ±0.01) cm, B area 
was (0.62 ± 0.07) cm, and C area was (0.48 ± 0.06) cm. The values obtained still meet the standards set by 
American Association of Physicist in Medicine (AAPM). 
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I. PENDAHULUAN 

Ultrasonografi (USG) termasuk dalam jenis pelayanan 
imaging diagnostic yang bersifat non-invasif sehingga tidak 
menyakitkan. Selain itu, pemeriksaan USG tersedia secara 
luas, mudah digunakan, lebih murah daripada metode penci- 
traan lainnya, dan relatif aman karena tidak menggunakan ra¬ 
diasi pengion [1]. Dengan beberapa kelebihan tersebut ten- 
tuny a kualitas citra menjadi faktor yang penting dalam hal pe- 
manfaatan USG, oleh karena itu perlunya dilakukan uji kon¬ 
trol kualitas dari USG, salah satunya terkait resolusi spasial. 


Resolusi spasial mengacu pada kemampuan melihat detail 
terkecil. Secara sederhana, resolusi spasial berhubungan de¬ 
ngan seberapa kecil sebuah objek dapat dilihat pada sistem 
pencitraan tertentu [2]. Detail kecil pada citra USG akan 
sulit sekali dianalisis menggunakan mata manusia karena be- 
rada pada orde pixel (elemen terkecil gambar), sehingga de¬ 
ngan keterbatasan itu diperlukan suatu metode yang dapat 
dijadikan sebagai instrumen menganalisis resolusi spasial. 
Belum adanya penelitian sebelumnya yang secara khusus di¬ 
lakukan untuk membuat suatu metode terkait proses analisis 
dari resolusi spasial hasil citra USG. 
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(a) Linear array 


(b) Curved array 


(c) Annular array 


Gambar 1: Tipe tranduser berdasarkan jenis array (diambil dari http://www.examrefresh.com/transducer-arrays/). 


Citra pada USG bergantung pada tranduser yang digu- 
nakan. Berdasarkan jenis tranducer arrays secara garis be- 
sar dapat dibedakan menjadi tiga, yaitu linear array, curved 
array , dan annular array [3], seperti ditunjukkan Gambar 1. 

Pada penelitian ini dibahas transduser dengan jenis curved 
array , dimana pada tranduser jenis ini akan menghasilkan 
citra yang memiliki bentuk seperti Gambar 1(b). Sistem koor- 
dinat pixel gambar objek yang dihasilkan dari transduser de¬ 
ngan jenis ini tidak dalam sistem koordinat kartesian karena 
menyesuaikan dengan permukaan dari transduser. Hal ini 
berdampak pada proses analisis resolusi spasial yang tidak 
dapat dilakukan dengan proses perhitungan secara sederhana, 
dengan kata lain perlu dilakukan transformasi citra tersebut 
menjadi koordinat kartesian agar lebih mudah dalam men- 
ganalisis resolusi spasial. 

Untuk melakukan uji kontrol kualitas diperlukan instrumen 
yaitu phantom ultrasound. Phantom ultrasound pada umum- 
nya memiliki harga yang cukup mahal sehingga tidak banyak 
tersedia untuk keperluan kegiatan belajar. Oleh karena itu 
diperlukan home made ultrasound phantom yang terjangkau, 
tahan lama, dan dapat digunakan lagi. 

Beberapa penelitian sebelumnya telah dibuat phantom ul¬ 
trasound dengan berbahan dasar gelatin, agar, gel parafin dan 
silicon rubber. Pada phantom berbahan dasar agar dan gelatin 
tidak dapat bertahan lama sehingga tidak dapat digunakan un¬ 
tuk jangka panjang, selain itu pada phantom dengan bahan 
dasar gelatin membutuhkan refrigerator dalam proses pem- 
buatannya [4] dan phantom dengan berbahan dasar gel parafin 
sulit dalam mengubah kecepatan suara dari material yang di¬ 
gunakan [5]. Pada phantom dengan berbahan dasar sillicon 
rubber memiliki keuntungan lebih tahan lama, dapat digu¬ 
nakan lagi dan aman digunakan [6]. 


II. METODOLOGI 

Penelitian ini menggunakan pesawat USG 3D merk Min- 
dray model: DP-10 dengan transduser model: 35C50EB. In¬ 
strumen yang digunakan untuk melakukan analisis resolusi 
spasial adalah phantom sederhana yang terbuat dari campuran 


8 cm 



Gambar 2: Desain phantom berisi objek berupa susunan paku 
berdiameter 2 mm dengan jarak antar paku di daerah (A) 9 
mm, di daerah (B) 6 mm, di daerah (C) 4,5 mm. 


sillicon rubber dan katalis. Proses analisis lebih lanjut meng¬ 
gunakan komputer dengan software pengolah citra. Desain 
phantom ditunjukkan pada Gambar 2. 

Tahapan untuk menganalisis resolusi spasial dimulai de¬ 
ngan dilakukannya scanning phantom dengan adanya pe- 
ngaturan depth dan gain. Depth merupakan salah satu tools 
pada mesin USG yang berfungsi mengubah bidang gambar 
yang ditampilkan dalam peningkatan gradasi satu sentime- 
ter. Menambah kedalaman akan mengurangi resolusi gambar. 
Pada kedalaman yang lebih tinggi, transduser perlu menem- 
puh jarak yang lebih jauh, oleh karena itu laju bingkai dan 
resolusi keduanya lebih rendah, sedangkan gain merupakan 
tools yang digunakan untuk menyesuaikan kecerahan keselu- 
ruhan dari gambar ultrasound real-time (B-mode). Perubahan 
pada gain dapat membuat gambar yang terlalu gelap menjadi 
lebih terang atau sebaliknya, sehingga dapat dilihat perbedaan 
tekstur yang halus. 

Selanjutnya dilakukan image processing dari citra hasil 
USG yang di dapat. Citra hasil USG dibaca menggunakan 
komputer yang telah ter-install software pengolah citra di¬ 
mana format citra dari RGB diubah menjadi grayscale. Tahap 
terakhir yaitu perhitungan resolusi spasial dan jarak antar ob¬ 
jek dilakukan dengan pendekatan persamaan lingkaran. 
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Gambar 3: Hasil citra USG untuk daerah B. 


ning, citra objek yang dihasilkan menunjukkan susunan objek 
berada pada lintasan yang melengkung. Oleh karena itu, 
perlu digunakan pendekatan persamaan lingkaran untuk dapat 
menghitung resolusi spasial arah tangensial dari citra USG 
tersebut. 


Perlu diketahui terlebih dahulu dimana letak pusat 
lingkaran menggunakan pendekatan persamaan lingkaran. 
Ditunjukkan pada Gambar 5, bahwa bentuk dari hasil citra 
USG menyerupai potongan dari sebuah lingkaran. Dengan 
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Gambar 4: Sistem koordinat pixel. 


III. HASIL DAN DISKUSI 

Penelitian ini awalnya dilakukan dengan terlebih dahulu 
membuat phantom sederhana yang terbuat dari bahan silicon 
rubber. Beberapa sample phantom telah dibuat dengan vari- 
asi konsentrasi katalis. Objek dapat nampak saat dilakukan 
scanning seperti pada Gambar 3, digunakan campuran katalis 
1% dari total campuran volume. Phantom yang digunakan 
sebagai instrumen untuk mendapatkan data berasal dari 
campuran 99 ml sillicon rubber dan 1 ml katalis. 

Tahapan dalam analisis dimulai dengan processing image , 
dimana gambar asli dari hasil citra USG yang berformat 
RGB diubah menjadi grayscale. Hal ini bertujuan untuk 
membuat warna gambar berada pada rentang hitam dan putih 
atau dalam aras keabuan disebut sebagai biner, agar lebih 
memudahkan untuk dapat membedakan antara background 
dan objek yang diamati. 

Sebuah gambar didefinisikan dalam fungsi 2-dimensi f(x, 
y), dengan x dan y adalah koordinat spasial dan amplitudo 
f di (x, y) tertentu, atau biasa disebut graylevel. Dalam 
pemrosesan gambar, gambar ini diwakili oleh pixel (elemen 
gambar) matriks M x N dengan sistem koordinat ditunjukkan 
pada Gambar 4 [1]. 

Proses analisis resolusi spasial dilakukan salah satunya 
dengan mengukur jarak antar objek yang berdekatan. 
Berdasarkan desain phantom yang dibuat, susunan antar 
objek berada pada kedalaman yang sama dengan perbedaan 
jarak berada pada arah horisontal. Setelah dilakukan scan- 


„ • 

^(xl,yl) 

• 

<*2,i2) 

• 

(*3,y3) 



4« 



Gambar 5: Pendekatan persamaan lingkaran pada citra USG. 


menggunakan beberapa titik acuan, pusat lingkaran tersebut 
dapat diketahui. 

Persamaan lingkaran dengan pusat P(a,b) dan jari-jari(r): 

r 2 = (x — a) 2 + (y- b ) 2 (1) 

Perhitungan nilai koordinat dari pusat lingkaran tersebut 
membutuhkan minimal tiga titik acuan dari gambar hasil 
citra yang memiliki jarak dengan pusat(jari-jari) sama besar. 
Pemilihan titik acuan ditunjukkan pada Gambar 4. Dari ketiga 
titik acuan tersebut akan didapatkan persamaan: 

r 2 = (xi - a) 2 + (yi - b) 2 (2) 

r 2 = ( x 2 - a) 2 + ( y 2 - b ) 2 (3) 

r 2 = (x 3 - a) 2 + ( y 3 - b) 2 (4) 


Dengan menggunakan metode substitusi akan dihasilkan persamaan baru untuk menemukan titik pusat (a,b): 


Qr 2 ) 2 - (*i) 2 + (y 2 ) 2 - (yi) 2 - 2bgn - vi) 

2(x 2 - x 2 ) 

[{x 2 ? - (*i) 2 + ( y 2 ) 2 - ( Vl ) 2 ] - ^[(x 3 ) 2 - (an) 2 + (y 3 ) 2 - (yi) 2 } 
2[{y 2 -yi)-{ x ^{y 3 -y 2 ))} 


( 6 ) 
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Setelah mendapatkan nilai titik pusat a dan b, dapat dihi- 
tung besar jari-jari (r) dari lingkaran tersebut dengan mema- 
sukkan nilai keduanya pada Pers.(2), (3), atau (4). Sebagai 
contoh, dalam Gambar 6(a) ditunjukkan citra dari sebaran tiga 
objek yang berhasil di-scanning dengan parameter frekuensi 
4,5 MHz, depth 3,3 cm, dan gain 168 dB. Pada Gambar 6(b) 
ditunjukkan pada arah radial dari tiap objek tersebut, terdapat 
tiga area yang dilalui. Daerah 1 menunjukkan area sebaran 
objek jika dilihat dari permukaan phantom , sedangkan daerah 
3 biru menunjukkan sebaran objek jika dilihat pada area yang 
lebih dalam dari area objek sebenarnya. Panjang busur K pada 
permukaan phantom diberikan oleh: 


K= - x 2 t tR 

360 


(7) 


Berdasarkan Pers.(7) panjang busur akan sebanding de¬ 
ngan jari-jari, sehingga semakin dalam areanya (semakin be¬ 
sar jarin-jariny a) akan terdapat penambahan jumlah pixel dari 
panjang busur. Meskipun jumlah pixel pada tiap-tiap panjang 
busur berbeda, akan tetapi dalam pengintrepretasian citra, 
daerah 1, 2, dan 3 harus menunjukkan panjang yang sama 
sebagai ukuran jarak objek pda kondisi nyata. Hal itu mem- 
buat proses analisis resolusi spasial arah tangensial citra USG 
tidak dapat langsung dilakukan dengan menghitung seberapa 
banyak jumlah pixel antara satu objek dengan objek yang 
berdekatan, perlu dilakukan proses transformasi. 

Terdapat satu variabel yang konstan meskipun terdapat 
perbedaan jumlah pixel untuk masing-masing panjang busur 
pada jari-jari tertentu, yaitu sudut yang dibentuk antara 
masing-masing objek dengan garis normal. Nilai sudut ini 
dapat digunakan sebagai pendekatan dalam menentukan nilai 
acuan perhitungan jarak antar objek. Pada citra USG sudah 
terdapat skala yang menunjukkan panjang area tertentu pada 
citra yang mewakili panjang pada kondisi nyata (dalam cen¬ 
timeter). Sebagai contoh, dalam Gambar 7 ditunjukkan area 
sepanjang panah kuning mewakili panjang 2 cm pada ukuran 
sebenarnya. Diasumsikan panjang dari skala itu sama dengan 
panjang pada busur di permukaan phantom, karena pada per¬ 
mukaan phantom yang dimungkinkan memiliki jumlah pixel 
yang sama dengan panjang skala yang diketahui. 

Pendekatan sudut untuk mendapatkan besar acuan (k) de- 
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Gambar 7: Skala pada hasil citra USG. 


ngan skala acuan (dalam cm) dihitung dengan persamaan: 

k = — ( 8 ) 

r 

dengan r adalah radius ke permukaan (dalam pixel) dan dida- 
patkan dari informasi acuan yang terdapat pada citra USG 
seperti yang ditunjukkan Gambar 8. Nilai dihitung dengan 
menggunakan persamaan: 

SA = y al - y a2 (9) 

Perhitungan jarak antar objek diakukan dengan scanning 
area di arah tangensial radiasi, yaitu dengan menentukan ter- 
lebih dahulu posisi radial objek dan sudutnya dengan per¬ 
samaan: 


Rib = (Xob - a) 2 + (Vob ~ b) 2 ( 10 ) 

, —in x ob ~ a 

tan 0 ob = - - (11) 

Vob ~ b 

kemudian dilakukan proses scanning objek yang tersebar sep¬ 
anjang area busur O. Untuk menghitung posisi objek sebe- 
namya (S ob ) dalam cm, digunakan persamaan: 

S 2 ob = 6 -f ( 12 ) 
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Gambar 8: Area scanning objek. 


TABEL I: Nilai hasil perhitungan. 
Daerah objek S 0 b(cm) JS*(cm) a oSE 


A 

0,85 

0,90 

0,30 

0,01 

B 

0,62 

0,60 

0,25 

0,07 

C 

0,48 

0,45 

0,15 

0,06 


*Jarak antar objek sebenarnya pada phantom 

kemudian didapatkan grafik yang menunjukkan hubungan an- 
tara intensitas dan posisi objek, seperti ditunjukkan pada 
Gambar 9, bahwa jarak antar puncak adalah jarak antar ob¬ 
jek sebenarnya. 

Berdasarkan grafik yang dihasilkan, dapat ditentukan dua 
nilai. Pertama nilai jarak antar obyek dengan menghitung 
jarak antar puncak tertinggi dan yang kedua nilai standar de- 
viasi a, dimana nilainya sebanding dengan setengah dari lebar 
spektrum (£). 


. - y 

(13) 

SE = 4= 

(14) 

\/n 


Pada daerah A, B, dan C dilakukan scanning area se- 
banyak 10 kali, sehingga n = 10. Lebar spektrum yang di¬ 
hasilkan cukup besar dikarenakan diameter dari obyek sebe- 
namya (dalam hal ini sebagai lebar obyek) berada pada orde 
mm, dimana orde ini sama dengan orde jarak antar obyek. 

Hasil dari perhitungan jarak antar objek untuk tiap daerah 
berdasarkan grafik yang diperoleh ditunjukkan Tabel I. 
Berdasarkan Tabel I, untuk masing-masing daerah dilakukan 
10 kali scanning , didapatkan nilai standar deviasi untuk 
daerah A sebesar 0,30, daerah B sebesar 0,25, dan daerah C 
sebesar 0,15. Untuk daerah A memiliki standar error sebe¬ 
sar 0,01, daerah B sebesar 0,07, dan daerah C sebesar 0,06. 
Berdasarkan standar error yang didapatkan, perhitungan jarak 
antar objek untuk daerah A sebesar (0,85 dk 0,01) cm, daerah 
B sebesar (0,62 db 0,07) cm, dan daerah C sebesar (0,48± 
0,06) cm. Dengan hasil perhitungan tersebut. 




Posisi, x(cm) 
(b) Daerah B 



Posisi, x(cm) 

(c) Daerah C 


Gambar 9: Grafik jarak antar objek. 


IV. SIMPULAN 


Berdasarkan hasil yang didapat, campuran sillicon rub¬ 
ber 99 ml dan katalis 1ml dapat digunakan sebagai phantom 
sederhana untuk melakukan uji kontrol kualitas USG dalam 
hal resolusi spasial dengan menghitung jarak antar objek yang 
berada pada arah tangensial dari radiasi menggantikan phan¬ 
tom buatan pabrik. Pendekatan persamaan lingkaran dapat 
digunakan untuk melakukan transformasi hasil citra USG 
melalui perhitungan sebaran obyek pada arah tangensial radi¬ 
asi dengan hasil perhitungan daerah A sebesar (0,85 ± 0,01) 
cm dari posisi sebenarnya 0,90 cm, daerah B sebesar (0,62 d= 
0,07) cm dari posisi sebenarnya 0,60 cm, dan daerah C sebesar 
(0,48± 0,06) cm dari posisi sebenarnya 0,45 cm. 
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